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Abstract 

Several different vertex detector designs for the Linear Collider Detector (LCD) are evaluated 
in the context of measurements of the branching ratio of the standard model Higgs particle 
decaying into charm/anti-charm final states. Fast Monte Carlo simulations are used to model 
the detector and neural network-optimized flavor tagging is used to perform the measurements. 
These tools arc used to study the effects of pixel resolution, material thickness, and inner 
layer radius on the flavor tagging efficiency-purity curve and ultimately on the branching ratio 
measurement error. 

1 Introduction 

A worldwide effort is now underway to develop a future Linear Collider and a corresponding Linear 
Collider Detector (LCD). Compromises may have to be made in the design of such a detector 
between its desired performance and other considerations such as the restrictions imposed by beam- 
related backgrounds in the interaction region as well as the total cost of the detector. In this paper, 
we use Monte Carlo simulations to model the performance of various detector configurations in 
the contexts of flavor tagging and sensitivity in the measurement of the branching ratio of a 120 
GeV Higgs boson to charm/anti-charm final states. Such a measurement would help determine 
whether the condensation of a single Higgs boson is sufficient to generate the masses of particles. 
The dependence of the measurement precision on the various detector parameters used provides 
useful data for optimizing the final design of the detector. 

The branching ratios BR(H° -> cc) and BR(H° -> bb) may be simultaneously measured using 
the SLD technique |Q] for measuring the the analogous quantities Rf, and R c in Z° decays. The 
branching ratios may be measured by performing flavor tagging and then counting the number of 
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6-tagged, c-tagged, and untagged decays. These data can then be unfolded using the flavor tagging 
efficiency matrix. 

Efficient and pure flavor tagging requires precise impact parameter resolution in order to dis- 
tinguish heavy flavor decay tracks from fragmentation tracks from the primary interaction point 
(IP). This tracking resolution may be written as ctq © &ms/(p ' sin 3//2 0) for a cylindrically sym- 
metric detector where 9 is measured from the z axis and p is the particle's momentum, do is the 
asymptotic resolution at infinite momentum and it scales linearly with the detector's single hit 
resolution. For finite momentum, the total resolution is usually dominated by the second term 
where &ms includes the effects of multiple Coulomb scattering and thus scales as the square root 
of the detector layer thickness. Both errors scale roughly linearly with the inner radius of the 
detector which is effectively the extrapolation distance from the innermost hit inward to the IP. 
The effects of variations of these three detector parameters on flavor tagging are described below. 
Most of the same assumptions are made as in a previous study of the same topic by the University 
of Oregon group [||, and our results are very consistent with theirs. However, the current study 
explores several additional detector configurations. 



2 The Analysis Tools 

The various vertex detector configurations are simulated using the LCDTRK programQ which, 
given an input detector geometry, generates smearing tables for the track parameter error matrices. 
Non-Gaussian measurement errors are not modelled. These smearing tables are input into the 
LCDROOT[|| package for the event generation and fast MC detector simulation. 

To calibrate the fast MC, the SLD tracking resolution is simulated by replacing the tracking 
chambers of the small detector (SDMAR01)[||] with a model of the SLD tracking chambers || 
0] [|| . This model is tuned to reproduce the measured SLD track position resolution of ~ (9 © 
33/((p/GeV) -sin# 3 / 2 ) fim. For this tuning, the beampipe thickness is increased from 0.52% X° to 
0.65% X° in order to accomodate the increased measurement error from Rutherford scattering tails. 
The magnetic field for the drift chamber is set to 0.6 T and interaction point (IP) measurement 
errors (4 • 4 • 20)/im are modelled by smearing the MC IP position. 20K Z° pole hadronic decay 
events are generated and the LCDJetFinder and LCDVToplGhost packages in LCDROOT are used 
to find jets using the Durham algorithm, and to perform topological vertexing j|] |uj within the 
jets. 



A jet flavor tag analogous to the SLD method] 11 1 is then developed. The tag proceeds in 
three stages, each of which is optimized with the Stuttgart neural network trainer. First, the best 
candidate b vertex is selected from the vertices found in the jet. This selection is based on the 
vertex flight distance, the flight distance normalized by the vertex measurement error, the angle of 
the vertex flight direction with respect to the jet axis, the number of tracks in the vertex, and the 
raw vertex mass computed from the track momenta. 

Next, tracks are attached to the b vertex using the following criteria: the transverse distance 
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of the track to the vertex flight axis at the point of closest approach (POCA), the longitudinal 
distance of the POCA on the vertex flight axis to the IP, this longitudinal distance divided by the 
vertex flight distance, the angle of the track momentum with respect to the vertex flight direction, 
and the track impact parameter normalized by its measurement error. The vertex kinematics are 
then recalculated using all of the attached tracks. 

Finally, the jet flavor is identified using the .Pr-corrected vertex mass, the vertex momentum, 
the vertex flight distance, the flight distance normalized by the vertex measurement error, the 
number of vertices found, and the number of I-prong vertices found. The results of the jet flavor 
tag are shown in figure |l|. These efficiency-purity curves are very similar to the the flavor tag 
results of the full SLD MC shown in figure ^. The fast MC b tag has slightly lower purity perhaps 
because the 'ghost track' vertexing algorithm used in the fast MC tag instead of the 'resolubility' 
vertexing algorithm used in the SLD tag. The 'ghost' algorithm is more susceptible to creating 
vertices from poorly measured fragmentation tracks, thus incorrectly increasing the vertex mass. 
The discrepancy is very small however, and we conclude that the fast MC using LCDTRK and 
LCDROOT is sufficiently good for our purposes. 

3 Flavor Tagging in Higgstrahlung Events 

To measure the various Higgs branching fractions, we look at Higgstrahlung events in which the 
final state Z° is identified by reconstructing its decay into e + e~ or fi + pT and a 120 GeV Higgs 
is identified through its recoil mass. The efficiency of the recoil mass technique to tag a 120 GeV 
Higgs is estimated to be ~ 31% with ~ 44% impurity from Z^Z® events due to degradation of the 
beam energy constraint from beamstrahlung and initial state radiation. To reduce the Z a — > cc 
background, a cut of | cos 9\ < 0.6 is imposed on the jet axes to remove a large fraction of the 
polarized Z°'s produced by the 80% polarized electron beam. The resulting Higgs efficiency is 
~ 19% with 75% purity. Since only leptonic decays of the Z° are used, the Higgs decay jets can 
be selected unambiguously in this study, but only ~ 6% of the Higgstrahlung events can then be 
used. However, jets from Higgs and hadronic Z° decays can be separated some fraction of the time, 
and so the fraction of usable events in a full analysis can very likely be increased by a factor of 4. 
Assuming a 60 fb Higgstrahlung cross-section, the total number of usable tagged Higgs decays in a 
500 fb -1 sample is then estimated to be: 

N H « 500 fb' 1 • 60 fb • 0.19 • 0.06 • 4 ps 1370. (1) 

For each detector configuration, 20K Higgstrahlung events are generated in which the Z is forced 
to decay into e + e~ or and the 120 GeV Higgs decays hadronically. Since the branching ratio 

BR(H — > cc) is expected to be low, to get a good statistical sample of charm jets, an additional 
20K Higgstrahlung events are generated in which the Higgs are now forced to decay to cc. The 
magnetic field is set to 5 Tesla. 
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The IP position measurement is modelled to have (4 x 4 x 6)fim resolution. These numbers 
are motived by the SLC IP position measurement []l2| which is composed of measurements of the 
beam position averaged over many events, and an event-by-event measurement of the longitudinal 
position along the beam of the primary vertex. The SLC beam position is measured by putting 
fragmentation tracks from 30 consecutive hadronic events into a common 2D vertex in order to 
average over the jet axis distribution in eft. At luminosities of 10 34 cm -2 sec -1 the event rate 
for W + W~ and Z ^/ events at the Linear Collider should be comparable to the Z° event rate 
(~ 200/hour) at the SLC, and this technique should still work provided that the beam does not 
move by more than a few microns per hour. The SLC longitudinal IP position is measured event- 
by-event using the point of closest approach to the beam of the median track in the track-to-beam 
impact parameter distribution of tracks in the event. The resolution of this technique is expected 
to be of the same order as the average track position resolution. The changes in the longitudinal IP 
resolution when using different detector configurations have been neglected and so the changes in 
the branching ratio measurement precision for the most extreme detector configurations in table || 
may be slightly underestimated. 

The vertex selection and track attachment neural networks are trained only once for the standard 
SDMAR01 configuration HI, but the jet flavor tagging neural net is retrained using these 40K 
events for each new detector configuration. 20K of the events are used for training and 20K for 
validation to avoid overtraining. The same training regimen is used for each configuration: 1000 
cycles of standard backpropagation with the learning rate set to 0.2 and then 1000 cycles each 
of backpropagation with momentum 0.5 and the learning rate set to 0.1, 0.05, and 0.01. Studies 
showed that additional training beyond this basic regimen yield only marginal improvements in the 
flavor tagging. Small improvements could also be made using a 300K event sample for the training, 
but working with such large sample sizes was not feasible for this quick study. 

Since each hadronic Higgs decay produces both a quark jet and an anti-quark jet, there are two 
chances to tag the decay type. The event level flavor tag may be defined as follows. If either or 
both jets are 6-tagged, then the event is 6-tagged. If either or both jets are c-tagged and neither 
is 6-tagged, then the event is c-tagged. Otherwise, the event is untagged and considered to be a 
decay into light quarks. For example, the event flavor tagging efficiencies E = Eij for identifying 
flavor i as flavor j may be expressed in terms of the jet flavor tagging efficiencies as: 

E C c — ^cc ^^cc ' £cu 

= 2-e cc -(l-e cb )-e 2 cc (2) 

Ebc = e bc + ^ e bc • tbu 

= 2 • £ bb ■ (1 - e bb ) - el c (3) 

Here, the unitarity of the jet flavor tagging efficiency matrix is used for the second equality and i 
and j take on values of b, c and u for untagged. Values of ebb ~ 0.6 and e c b ~ 0.02 are typical. In 
these formulae, the effects on flavor tagging due to energy correlations between the two jets in a 
decay are neglected. However, since the BR{H — > cc) measurement is dominated by the H —* bb 
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background, these correlations can only improve the measurement by reducing E^c by making it 
less likely that the two b jets are both tagged as c. 

4 Measuring the Branching Fractions 

The branching ratios BRh = (BR H ^ bb , BRh-+ C c, BRh^uu) m ay be extracted by counting the 
number of event level flavor tags of each flavor T = (T&, T c , T u ) by solving the matrix equation: 

T = E (BR H • N H + BR Z • N z ) (4) 

to get: 

ABR H = E" 1 • (T - E • BR Z • N z ) /N H (5) 

Here, Nh and Nz are the numbers of Higgs and Z 0, s in the Higgs-tagged sample, and the Z° decay 
branching ratios BRz are defined in the same manner as BRh- The flavor tagging efficiency is 
assumed to be largely independent of whether a jet originated from a Higgs decay or from a Z° 
decay. This assumption may need to be revised when considering Higgs masses much larger than 
the Z° mass. 

Neglecting the error in the efficiency matrix (measured by counting the mixed tag rates where 
one jet is tagged as flavor i and the other as flavor j) and the error from subtracting the Z° 
background (presumably measured from the shape of the recoil mass distribution), the branching 
ratio error may be estimated as: 

ABR H « E- 1 • VT/Nh (6) 

In particular, if the Poisson fluctuations AT], and AT U are ignored then for Higgs to cc decays: 

ABRs^ce « E- 1 • ^/T C /N H (7) 

The number of c-tagged events is estimated as 

T c w (E cc ■ BRh^cc + E bc ■ BR H ^ b i) ■ Nh + E cc ■ BRz^cc • Nz (8) 

where the small contamination from BR z ^ bb is ignored. The fractional measurement error is then: 

ABRh^cc ^ \J E cc ■ BRh^cc + Eb c ■ BR H ^ bb + E cc ■ BRz^ c s • Nz/Nh . , 

BRh^cc E cc ■ BR H ^ c5 ■ y/IVF 

It is estimated that the Higgs tag has 75% purity and so N z /Nh « 0.33%. It is also assumed that 
the true values of the branching ratios are: BRh^cc = 3%, BR H ^ bb = 72%, and BRz-> C c = 12.4%. 
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In figures |5|- 1| below, for various vertex detector configurations listed in table [l], the efficiency- 
purity curve of the event flavor tag is plotted. These curves are generated by varying the value of 
the cut on the neural network output. The scaled precision P, defined as: 

P = 0Va" • ABRh^cc/BRh^cc (10) 

is also plotted as a function of the event flavor tagging point chosen on the efficiency-purity curve. 
This quantity is a measure of the analyzing power and can be divided by y/Nn ~ v 7 !^) ps 37 for 
a 500 fb _1 data set to get the fractional branching ratio error. 

For example, for the standard SDMAR01 configuration with 5 fim hit resolution and 0.12% X° 
thickness and inner radius 1.2 cm, figure || indicates that the best precision occurs at E cc ~ 0.66 and 
Ejj C ~ 0.06 and the relative branching fraction error can be calculated to be ABRh^ c£ / BRfj^ c5 ~ 
13.8/37 ~ 0.37. For Nh ~ 1370 and Nz ~ 452, the breakdown of the tagged charm events is ~ 27 
correct tags from H° — > cc, ~ 59 mistags from H° — > bb, and ~ 36 correct background tags from 
Z° — > cc. The mistags from Z° — > bb are neglected. Our result is very consistent with that of the 
Oregon result M of ABRh~^ C c/ BRh^cc ~ 0.39 for the same model. 

5 The Results 

The single hit resolution a res is varied (figure ||) between 1/iin and 10/rni to cover resolutions 
ranging from CCD pixel detectors and of silicon strip detectors. For example, changing from the 
5//m resolution SDMAR01 detector to a 1/iin detector only decreases the error by a factor of 
13.3/13.8 = 0.964. The changes in precision may be estimated as AP '/ Aa res ~ 0.12 fim . 

A large range of variations is chosen for the detector layer thickness t, from 0.03% A to 1.00% A 
again to model anything from the thinnest CCD structures to silicon strips (with readout electronics 
electronics on the ends). These variations, shown in figure |6| resulted again in changes of precision 
of magnitude similar to those of the single hit resolution variations. These changes may be roughly 
parameterized as AP/ Ay/i ~ 1.38 where t is measured in %A°. The larger effects of simultaneous 
variations in both single hit resolution and layer thickness are shown in figure 

The dependence on the vertex detector inner radius Rq is shown in figure ||. Not much difference 
can be seen, and the dependence is approximately AP/ARq ~ 0.17/cm. This result indicates that 
it may be possible to increase the detector inner radius to better avoid the e + e~ pair backgrounds 
without sacrificing much precision in the branching ratio measurement. 

Finally, the dependence of the measurement on the center-of-mass energy is shown in figure ^ 
for y/s = 250, 500 GeV for two different detector configurations. Changing the energy appears to 
have only small effects on the flavor tagging contribution to the measurement precision. Perhaps 
more important are the effects of different Standard Model backgrounds to the Higgs tag at the 
two different energies which have not been modelled in this study. 
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Table 1: The different detector configurations based on variations of SDMAR01 and the cor- 
responding track position resolution, do © ctms/(P ■ sin 3//2 6). The last entry represents the SLD 
model which has an enhanced 0.65% X° beampipe at 2.5 cm radius and has the SDMAR01 tracking 
chambers replaced with a model of the SLD vertex detector and central drift chamber. All other 
models have a 0.04% X° beampipe at 1.0 cm radius. 

6 Conclusion 

In this study, we have found that the precision of a BR(H° — ► cc) measurement is only modestly 
dependent on the variations of several detector parameters. This small dependence can be explained 
by the fact that the base configuration (SDMAR01) already provides exceedingly good tracking 
resolution. Even the worst detector configuration modelled offers resolution much better than 
that of any existing vertex detector. The modelled tracking resolutions (shown in table |l|) are 
much smaller than the distance scales given by the lifetimes of weakly decaying heavy hadrons, for 
example ctd ~ 140/mi. Flavor tagging analyses are therefore able to cut deeply into the exponential 
impact parameter distribution to achieve a high efficiency for separating heavy flavor decay tracks 
from fragmentation tracks. Further improvements in detector precision can only hope to recover 
the small remaining fraction of the exponential distribution at low proper time. 

It therefore does not seem necessary from the standpoint of flavor tagging to attempt to go to 
extraordinary lengths to achieve further improvements in resolution. However, other applications 
of precision tracking still need to be studied. For example, in measurements of parity violation 
or CP violation in new physics sectors or even at the Z° pole in a giga-Z run, it is necessary to 
determine not just the particle species, but also its charge. Even if a newly discovered particle does 
not have a weak decay scale lifetime, its charge may be measured if it decays to a b or c hadron. 
In either case, the loss of even a single secondary decay track creates an impurity in the charge 
measurement. It is therefore expected that these types of measurements may be more sensitive to 
changes in the tracking resolution. 
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Figure 1: The efficiency-purity curves of the LCDROOT fast MC jet 6-tag (left) and c-tag (right) 
for the SLD detector. The 'ghost track' algorithm is used for vertexing. 




c efficiency b efficiency 

Figure 2: The efficiency-purity curves for the full MC jet c-tag (left) and 6-tag (right) for the SLD 
detector. The 'resolubility' algorithm is used for vertexing. 
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Figure 3: The b — > c mistag rate versus 
the c — > c correct tag rate of the jet flavor 
tag for different numbers of neural network 
training cycles. 
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Figure 4: The efficiency-purity curve of 
the jet flavor tag for neural networks trained 
with 300K Z° events and with 40K Z° pole 
events. 
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Figure 5: (Left) The event mistag rate E^ c versus the event correct tag rate E cc for vertex detectors 
of different hit resolution but the same layer thickness of 0.12% X° and the same inner radius of 
1.2 cm. (Right) The corresponding fractional BR(H — > cc) precision. 
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Figure 6: (Left) The event mistag rate E\, c versus the event correct tag rate E cc for vertex detectors 
of different layer thickness but the same hit resolution of 5 fim and the same inner radius of 1.2 
cm. (Right) The corresponding fractional BR(H — ► cc) precision. 
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BR(H->cc) sensitivity, ZH(120)@500 GeV | 
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Figure 7: (Left) The event mistag rate E^ c versus the event correct tag rate E cc for vertex detectors 
of different layer thickness and different hit resolution but the same inner radius of 1.2 cm. (Right) 
The corresponding fractional BR(H — > cc) precision. 
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Figure 8: (Left) The event mistag rate Ef, c versus the event correct tag rate E cc for vertex detectors 
with different inner radius, using the same hit resolution of 5 /im and the same layer thickness of 
0.12% X°. (Right) The corresponding fractional BR(H — » cc) precision. 
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BR(H->cc) sensitivity, ZH(120)@250,500 GeV 
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Figure 9: (Left) The event mistag rate Et> c versus the event correct tag rate E cc at two different yfs 
for two detector models. The inner radius is kept at 1.2 cm. (Right) The corresponding fractional 
BR(H — > cc) precision. 
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Charm jet tagging, ZH(120)@500 GeV 
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